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Abstract. We study the number of periodic solutions in two first order non- 
autonomous differential equations both of which have been used to describe, among 
other things, the mean magnetization of an Ising magnet in the time- varying external 
magnetic field. When the strength of the external field is varied, the set of periodic 
solutions undergoes a bifurcation in both equations. We prove that despite profound 
similarities between the equations, the character of the bifurcation can be very different. 
This results in a different number of coexisting stable periodic solutions in the vicinity 
of the bifurcation. As a consequence, in one of the models, the Suzuki-Kubo equation, 
one can effect a discontinuous change in magnetization by adiabatically varying the 
strength of the magnetic field. 



PACS number: 05.50.+q 
1. Introduction 

In many lattice models the mean field approximation leads to ODEs with periodically 
time dependent right-hand side. For example, the Curie- Weiss model with Glauber spin 
dynamics in the thermodynamic limit leads to the so-called Suzuki-Kubo equation (see 
the derivation in [HI Ch. 3] or in fTJEj, and a numerical study in PQ). 
The Suzuki-Kubo equation is 

e-j- — — m + tanh[/3(m + ^,cos(27rr))], (1) 

GST 

where m is the average magnetization of a sample, e is the relaxation time, h is the 
amplitude of the applied magnetic field, and (3 is 1/kT, with T the absolute temperature. 
A seemingly similar first order non-autonomous equation, 

dec 

e— = ax + bx 3 + /icos(27rt), a, b el, (2) 

(I I 
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has been used as a simple generic model for switchable bistable systems. If a > 
and b < 0, it describes the overdamped dynamics of a particle in a quadratic potential 
with periodic forcing. Such equations have been studied in the context of laser optics 
(longitudinal mode instabilities in a semiconductor laser); see |H] for analysis and 
references to the optics literature. 

The issue we want to address is the number of 27r-periodic solutions for various 
values of h and other parameters. In particular, we will use h as the bifurcation 
parameter. This issue is of physical importance, since one wants to determine whether 
it is possible, say, in the Suzuki-Kubo context, to effect a discontinuous change in 
magnetization by continuously varying the amplitude of the applied magnetic field h 
(i.e. a first order phase transition). 

To compare the two equations effectively, we rescale the variables in (jlj in the 
following way: set, in sequence, 



x 



m + h cos(27rr), h = h\/l + 47r 2 e 2 , t = r H arccos ( — j . (3) 

27r Vvl + 47r 2 e 2 / 

Then in the new variables we have 
dx 

e— = — x + tanh/?x + h cos(27rf). (4) 

(JjL 

Comparing the two equations, we see that the nonlinearity in (Jlj) is sublinear while 
it is superlinear in (J2J). Furthermore, tanh(z) is not an entire function. Note also that 
the Taylor series expansion of — x + tanh f3x is 

B 3 

-x + tanh/3x= (/?- l)x-— x 3 + 0(x 5 ). (5) 

3 

Hence superficially it looks that the behaviour of solutions of (jlj) should be similar to 
that of (J2J) with a = (3 — 1 and b = —f3 3 . In fact, this is the implicit assumption 
under much of the physical literature. This assumption is incorrect. It turns out that 
for certain values of (3 and h equation (Jlj) has 5 periodic solution, a situation which is 
impossible for equation (J2J). The numerical simulations confirming this were reported 
in, for example, [3, Ch. 7]. Here we aim to establish this fact analytically. 

The structure of the paper is as follows. After collecting preliminaries below, in 
section 01 we perform the Liapunov-Schmidt reduction of the general equation, 
dx 

— = f(x) + hcos(2nt), (6) 

The reduction is complicated by the fact that we do not know the bifurcating solution 
explicitly. 

In section @] we use the results of the reduction to conclude that for small (3 
behaviour of (jlj) is similar to that of (J2J). Then we consider the case of large (3 in (0J) by 
explicitly examining the case of (3 = oo. The limiting equation is in fact a differential 
inclusion which we analyse explicitly, finding the range of h for which 5 solutions co- 
exist. Appealing to a continuation theorem j2], we conclude that such behaviour persist 
for large, but finite, values of (3 as well. 




Everywhere below we assume that the function /(■) is sufficiently smooth. As can 
be seen from (JHJ), we fix e = 1, as this does not affect any of the results. However in 
our numerics we take e of the order of 0.05 which "slows down" the time and results in 
nicer plots. 

Whenever we sketch the solutions of the equation we do it in (x, hcos(2irt)) plane, 

1. e. we plot the curve (x(t),y(t)) where 

x = f(x) - y, ^ 
y = h cos(27rt). 

This makes it easier to identify the periodic solutions (they are closed curves) and 
produces familiar hysteresis loop pictures for large h, figure ^ 

2. Background information 

The examples of function f(x) that we are going to consider here, f(x) = tanh/3x — x, 
(3 > 1 and f(x) = 72 — x 3 , 7 > share many important properties. Both are 
odd functions and the equation f(x) = has exactly three solutions. If we consider 
equation (JHJ) for h = 0, there will correspondingly be three stationary solutions. The 
central stationary solution x(t) = will be unstable and the other two are locally 
asymptotically stable. 

It follows from the implicit function theorem that for small h, there will be three 
periodic solutions, two stable and one unstable. Below we show that this situation 
persists at least until ho = max I> o f(x). Also we will show that for large h there is only 
one, stable, solution. The main aim of this paper is to uncover the bifurcation picture 
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Figure 2. Two possible minimal scenarios of the emergence of a single periodic 
solution. In scenario (a) three solutions merge in a subcritical pitchfork bifurcation. In 
scenario (b) the central solution first undergoes a supercritical pitchfork bifurcation, 
emitting two unstable solutions. These unstable solutions disappear in fold bifurcations 
upon meeting the stable solutions. 



for the intermediate values of h. The simplest scenario where three solutions can merge 
into one is a subcritical pitchfork as shown on figure |2fa). There, for simplicity, we 
sketch the means of periodic solutions, x = x(t)dt, versus h. 

Contrary to the above expectations, numerical investigations of the Suzuki-Kubo 
equation (jlj performed in [T3j (see also [SI Ch. 3] and j2]) show that for large values of f3 
and the parameter h, the equation has three stable periodic solutions (and therefore five 
periodic solutions altogether). The bifurcation diagram would then look like figure E[b). 
For small values of /3, JH| reported two stable periodic solutions at most, a situation 
compatible with figure a). The critical value of (3 at which the bifurcation picture 
changed was referred to as the tricritical points (TCP). 

The existence of TCP was disputed by [H] (see also P), who studied equation (J2J) 
numerically and reported that, although the stable non-symmetric solution changed fast 
in the vicinity of bifurcation, it was not disappearing in a "blue-sky catastrophe" as in 
figure I2jb). 

Our investigation concludes that, from mathematical point of view, both results 
are correct: there are at most 3 solutions in (J2J) and there is a TCP in equation (0J. 

2.1. Periodic solutions for small h 

Theorem 2.1. Let f(x) be a continuously differentiable function satisfying 



Put ho = min(/(a), —/(&)). Then for any h < h there is a unique periodic solution to 
equation (Q|) satisfying x(t) e [a, b] for allt. This solution is stable. 

Proof. The idea of the proof is to construct such a trapping region that any solution 
that enters it must stay there. Then we can appeal to a fixed point theorem to infer 




and 




(8) 
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that there is a periodic solution inside this region. Analysing stability of this solution 
we find that it must be stable, therefore it is unique. 

The function f(x) is monotone therefore there are unique solutions to the equations 
f(x a ) = ho and f(xf,) = —h . Since, for any t, x(t) = x a implies x'(t) = f(x a ) + 
/icos(27rt) > ho — h > 0, any solution to the right of x a will remain there. Similarly, 
solutions to the left of Xb will never increase past Xb- Thus [x a , Xb] is a trapping region 
and there is at least one periodic solution. 

Linearizing equation (0) shows that a solution xo(t) is stable if 

A = fo l % {£{s))dS (9) 
is negative. But since f'(x) < for all x e [x a , every periodic solution is stable. □ 

Corollary 2.2. Let x m be such that f(x m ) = max x> of(x), where f(x) is either 
tanh/Sx — x or jx — x 3 . Then for h < f(x m ) equation (GJ) has exactly three solutions, 
one in each of the intervals (— oo, — x m ), (—x m ,x m ) and (x m , oo). Two of these solutions 
are stable and one is unstable. 

2.2. Uniqueness of periodic solution for large h 

For large h we have the following general theorem, which, in particular, is valid for both 
functions f(x) that are of interest to us. 

Theorem 2.3. Let f{x) be a continuous almost everywhere differentiate function. If 
its derivative satisfies 

f'(x) < a < for \x\ > xo and f'( x ) < P f or oil x (10) 

for some a, (3 and xq, then, for sufficiently large h, equation has a unique stable 
periodic solution. 

Proof. Existence. Property implies that f(x) — > =Foo eventually monotonically as 
x — > ±oo. This means, in particular, that the equations f(x) = h and f(x) = —h have 
exactly one solution each for sufficiently large h. Denoting these solutions by Xh- and 
Xh+, we notice that the interval [xh-, Xh+] is trapping for the trajectories of JTUJ). Hence 
there must be at least one periodic solution. 

Uniqueness. We will prove the stability of all periodic solutions which will imply 
that there can be at most one. Let x(t) be a periodic solution, which is stable if A, 
defined by equation Q, is negative. Denoting the time spent by the solution inside 
[— xo,xo] by t + , we can estimate A by 

A < (3t + + a(l - t+). (11) 

We aim to show that the time t+ — > as h — ► oo and hence A < 0. 

To do it we define three regions in the plane (x,hcos(2irt)) which add up to the 
strip \x\ < x , see figure El The region A is defined by |/icos(27rt)| < y , the region B 
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Figure 3. Sketch of an orbit crossing [— xo,xo] and the three regions dividing the 
strip. 



by /icos(27rt) > y Q and the region B' by /icos(27rt) < — y . For now, we leave the choice 
of yo open. 

From the definition of region A, the time spent by the solution in A is at most 



t A = 2 i-arcsin(|). 



(12) 

2tt \hJ K J 

If M is the maximum of \f(x)\ on [— x ,x ], then \x'(t)\ > y — M and the time spent 
by the solution in B or B' during one visit is at most 

= (13) 
y - M 

It is easy to see that a periodic solution can visit B at most once during one period (and 
at most twice if it is not periodic); the same applies to region B' . Thus we conclude 
that the total time spent in the strip \x\ < Xq is 

Ax 



U < 



1 • {Vo\ 
— arcsm — 

7T V h / 



yo-M y ' 

Now if we take y to be equal, for example, to y/h, the right hand side will tend to zero 
as h — > oo. □ 
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3. The Liapunov-Schmidt reduction 

In this section we assume that function f(y) in equation (jHJ) is odd. It immediately 
follows that if x(t) is a solution of then so is x(t) = —x(t+ 1/2). In particular, there 
is always a "symmetric solution" which satisfies x(t) = —x(t + 1/2). 

We shall perform the Liapunov Schmidt reduction (LSR) for equation (JOJ) aiming 
to determine the character of the bifurcation of the symmetric solution as we change 
the parameter h. We shall denote by h cs the corresponding§ critical value of h. 

We start by defining the operator $ by 

=±-f(x) -hcos(2irt). (15) 

The periodic orbits are the zeros of this operator in the space of continuously 
differentiable 1-periodic functions C*([0, 1]) contained in L 2 ([0,1]). The reduction, 
based on a particular solution Xo(t), leads to the construction of a reduced function 
g: R x R — > R such that the solutions of g(x,h) = are locally in one-to-one 
correspondence with the solutions of (fIK| . with the solution Xo(t) corresponding to 
the zero solution of g(x, h) = 0. It is rarely possible to compute g(x, h) explicitly, but 
one can examine the bifurcation picture by computing the derivatives of g. 

First of all, however, we describe the symmetric solution xo(t) at the critical point 



3.1. The bifurcation condition 

The necessary bifurcation condition is that the differential of $ has non-zero kernel. 
This is because if Ker(o?$(x )) = {0}, the solution x of the equation $(x , hcs) — can 
be uniquely continued by the implicit function theorem. In our case the differential of 
$, which we denote by L, is 



d$(x)K) = <f>( x + £ Vl ) 



df 

= vx- j-{x)vi = Lvt (16) 

£=0 ° X 



Therefore, the necessary condition for bifurcation at xq = xo(h,t) is that there is a 
non-zero 1-periodic solution to 
df 

v = -g^( x o) v - ( 17 ) 
The general solution of (fTTj) is 

-t 

dx 



1 df 

v(t) = C exp I / — — (x (s))ds 



This solution is 1-periodic iff 

£^L( Xo (s))ds = 0, (19) 

which we will refer to as the bifurcation condition. 
§ the subscripts "cs" stand for "critical symmetric" 
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To proceed with the LSR we need to find a basis for KerL and (Range L)- 1 . The 
former is spanned by the function v(t) defined by (fTSj) with, for example, C — 1. The 
latter satisfies 

= (v*, Ly) = (-j f v* - |j>o(t)K, y) (20) 

for all y. Solving equation 

v* = ~%(x )v* (21) 
ox 

and again setting the arbitrary constant equal to one, we obtain 

v*(t) = exp (- £ ^( Xo (s))ds) . (22) 

Note that with our choice of arbitrary constants v(t)v*(t) = 1 which will be used often 
in what follows. 

3.2. Symmetries of the reduced function 

The function g(x, h) inherits the symmetries of the original operator $(x). In our case, 
the following proposition from [6, 4.3] holds 

Proposition 1. Let dim Ker L = dim RangeL = 1 and let KerL = span{v}. If there 
is a symmetry operator R: C^QO, 1]) — > Cp([0, 1]) which satisfies 

R 2 x = x, ®(Rx) = R®(x) and Rx (t) = x (t), (23) 

for all x G Cp([0, 1]) then Rv{t) is equal to either v{t) or —v(t). In the latter case the 
reduced function is odd in x: g(—x, h) = —g(x, h). 

It is easy to check that the symmetry operator satisfying conditions (J23j) is 
R<fi(t) = —4>{t + 1/2). Since v (t) is an exponential (see (JIBJ)), its sign does not change 
over [0,1] and Rv{t) = —v(t+ 1/2) must be equal to —v(t). Therefore our reduced 
function is indeed odd. 

As a consequence, we immediately get 

g(0,h) = Q, g xx (0,h) = 0, etc. (24) 

In addition, g x (0, h cs ) = since h cs is critical. Thus, to investigate the pitchfork 
bifurcation of the symmetric solution of (jUJ), if the bifurcation is indeed a pitchfork, 
it suffices to study ^(0, h cs ) and g xxx {0, h cs ). 

3. 3. Bifurcation function 

To study the remaining derivatives of g(x,h) we compute higher derivatives of the 
operator $ (here Vi,v 2 ,v 3 G Cp([0, 1])): 

d 2 ^(x)( Vl ,v 2 ) = -^l{x) Vl v 2 (25) 

d 3 f 

d 3 $(x)(v 1 ,v 2 ,v 3 ) = - -^(x)v 1 v 2 v 3 (26) 
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— $ = -cos(27rf) (27) 

oh 

=0. (28) 
We have already concluded that 

9k = ("*• ss*) " °- (29) 

g„=(v', <?*(»,»A = 0. (30) 



These formulae will come in handy below. 
The derivative ghx is given by 

»■ = (••■'' (I*) - <p *("- £ " lB (!*)))• (31) 

where i£ is the projection to Range L, 

Ey = y-(v\y)^— (32) 

If II 

In our case we have d (Jr$) = and, since (v* , J^3>) = = 0, the projector i? leaves 
invariant. Taking into account that v*v = 1 we get 



1 ^2^ 

^(x (s))m(s)c?s, (33) 



where u(s) = L 1 (Jr3>) is the solution of 
df 

u — — — (xo(s))u = — cos(27rt). (34) 
Ox 

Differentiating the identity <&(xo(h,t)) = with respect to h, we notice that u(s) = 
—dx (h,s)/dh. Eventually we get 

d f 1 df 

9hx(0,h cs ) = -— J —(x (s))ds, (35) 

therefore, at the stability-gaining bifurcation we should have 

g hx > 0. (36) 

Remark 1 . One consequence of theorem l2.3l is that the symmetric solution has to undergo 
at least one bifurcation in which it turns from being unstable to being stable. As of now, 
we are unable to prove that it undergoes exactly one bifurcation which is equivalent to 
showing that ghx is strictly greater than zero at all bifurcations. 

The character of the stability-gaining bifurcation is determined by the sign of g xxx , 

Theorem 3.1. If g XX x{®, h cs ) is positive (corresp. negative), the stability- gaining 
bifurcation from the zero-mean solution is a subcritical (corresp. supercritical) pitchfork. 
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Proof. We expand g(x, h) in Maclaurin series in x, taking into account that the function 
is odd, 

3 

g(x, h) = g x (0, h)x + g xxx (0, h)^- + 0(x 5 ) (37) 

o 

To understand the bifurcation picture, it suffices to take into account only the first two 
terms. If g xxx (0, h) > 0, the equation 

3 

g x (0,h)x + g xxx (0,h)^- = (38) 

o 

has 3 solutions if g x (0, h) < and 1 solution otherwise. Since g xxx (0, h cs ) > it remains 
so in a small neighbourhood of h = h cs . Thus the main question is the sign of g x (0, h). 

To prove that, as h crosses h cs , g x (0, h) goes from being negative to being positive, 
we can show that || 

9,(0, h ) = ~J %( x o(s))v{s) ds = -X (39) 

which is rather involved. Alternatively we can quote a result connecting the sign 
of g x (x,h) with the stability of the corresponding solution of the original differential 
equation [Bj. In our case, the solution corresponding to (0, h) is the symmetric solution 
which is stable (unstable) whenever g x (0,h) is positive (corresp. negative). Since the 
bifurcation is stability-gaining, we get three solutions for h < h cs and one solution for 
h > h cs . □ 

To study the sign of g xxx we use the following formula 

g xxx = (V , d 3 $(v, v, v) - 3rf 2 $ (v, L- l Ed 2 ${v, vj) ) (40) 

From the definition of the projector E, see ()32|1 . and identity (|30j) we conclude that 
Ed 2 $(v,v) = rf 2 $(t>,t>). Further, inverting L we get 

rt r tg^ 

dx 2 

First we evaluate the second part of (j40j). 
-s(v*, d 2 $ (v,L- x EdP$(y,v)) 

= " 3 J! V ^U {x)v2{t) S (X(S)MS) ds I dt ' (42) 
which can be written as 

3w 2 (t) 



L- 1 d 2 ${v,v)=v{t) I $(v,v)v*(s)ds = -v(t) I *p-{x{s))v(s)ds. (41) 



w'(t)w(t) dt 



i 



= 0, (43) 

J o 



2 

where w(t) = J * ^(x(s))v(s) ds. The integral above is equal to zero since w(l) = g. 
0. 

1 1 of which l|35|) is a consequence but not a proof 
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9xxx = (v*, d 3 ${v,v,v)} = - I ^{x {s))v 2 {s)ds. (44) 

Corollary 3.2. The stability- gaining bifurcation of the zero-mean solution of (0) with 
f(x) = x — x 3 is a subcritical pitchfork. 

Proof. In this case = —6 and, since v 2 is non-zero, g xxx > 0. □ 

A different proof of this (local) theorem for cubic / is given by Byatt-Smith [S]. 
That in the optical bistability equation (J2J) one can never have more than three 1-periodic 
solutions is proved very elegantly in [TU] . 

The above results will be used in section 14.11 to show that locally the situation in 
the Suzuki-Kubo equation (jlj) is similar to (J2J) if (3 is sufficiently small. 

4. The Suzuki-Kubo equation 

In this section we study the case 

f(x) = x - tanh(/3x) (45) 

in detail. 

We aim to demonstrate that for sufficiently large f3 there are values of h for which 
there exist at least 5 periodic solutions of (J3J. We also show that for small (3 the 
bifurcation of the zero- mean solution is a subcritical pitchfork (figure ED(a>)) which means 
that the maximum number of solutions is 3. The boundary between these two regimes 
is known as the tricritical point (TCP). 

4-1. The case of small ft 

Theorem 4.1. There exists (3q such that for all 1 < /3 < 0o and the zero-mean periodic 
solution xo(t) satisfying Mf^l we have 

J*^(x (s))v 2 (s)ds<0. (46) 

It is not hard to get the rough estimate, /3 = 4/3, which is sufficient for our 
purposes. 

Proof. Put tanh(/foo(s)) = C( s )- From fjlTj) . we have 

/ ^(x (s))v n (s) ds= [ v^^vis) ds = for all n > 1, (47) 
Jo dx Jo 

since v(t) is 1-periodic. Differentiating (J43j) . 

^L( Xo (s))={3-l-P( 2 ( S ), (48) 
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multiplying it by v 2 (s) and choosing n = 2 in (|47|) . we have 

1 C{s)v 2 {s)ds= ( ^- I v\s)ds. (49) 



o P Jo 

For the third derivative of f(x) we now get 

— (x (s))v 2 (s) ds 

= -2f3 3 f v 2 (s)ds + 8(3 3 f v 2 (sK 2 (s)ds-6f3 3 f (\s)v 2 (s)( ds 
Jo Jo Jo 

<-2f3 3 [ v 2 (s) ds + 8f3 3 [ ( 2 (s)v 2 (s)ds 
Jo Jo 

= 2/3 2 (3 / 3-4) f v 2 (s)ds, (50) 
Jo 

where we discarded the integral of a non-negative function and then used (|49jl . Thus 

the integral on the left-hand side of (f5T?j) is negative if (3 < 4/3. □ 

Remark 2. It is easy to improve the above estimate to 3/2 if, instead of discarding the 
integral of ( A v 2 altogether, we deduce from (f4*%j) 

4-2. The P = oo case 

In this section all assertions are verified by explicit computation. 

When we take (3 — oo, instead of the ODE (j3J) we are left with the differential 
inclusion 

dx 

e— G —x + Sgn(x) + h cos(27rt), (52) 
dt 

where Sgn(x) is the multivalued sign function, 

-1, if x < 0, 

Sgn(x) = { [-1,1], ifx = 0, (53) 
1, if x = l. 

Note that if a periodic solution x(t) has no internal zeroes, it solves a simple linear 
ODE and can be computed explicitly. For example, a for h sufficiently small, a positive 
periodic solution is given by 

., „ , cos(27rt) + 2ire sin(27rt) 
^ = 1 + h 1 + 4^ ' (54) 

from which we conclude that this (asymmetric) solution touches zero when 



h = h ca = Vl +4vr 2 e 2 . (55) 
For example, if e = 0.05, h ca = 1.048187027 to 10 d.p. 




Figure 4. Evolution of the unstable zero-mean solution of inclusion (|52l) for h < h cs . 
Top left: solution is identically zero when h < 1. Top right and bottom left: emergence 
and growth of nontrivial parts. The solution will experience a bifurcation when the 
trivial parts, where x(t) = on an open interval, disappear completely. This happns 
at h = h cs (bottom right). 



Next we want to understand the bifurcation^ of the zero-mean solution, figure |U 
When h < 1, the zero- mean solution is identically zero and is unstable. After h crosses 
1, the solution gets two symmetric non-trivial parts, but there are still two intervals in t 
on which x(t) is identically zero. The bifurcation happens when each of these intervals 
collapses to a point. 

Let h cs be the unique value of h for which there exists a periodic solution x c (t) of 

a complete bifurcation picture and animated periodic solutions described in this section are also 
available at http://www.math.tamu.edu/~berko/bistable 
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A hcos(27it) 



h cos(2tu t) 




Figure 5. The evolution of the right unstable solution until its annihilation at h = h c 



the differential inclusion with the property that there is to such that x' c (to) = 1 and 
x' c (to + 1/2) = —1, see figure ^bottom right). 

An involved but elementary computation (see Appendix A ) shows that this happens 

for 

,_, ^ y/g + 47r 2 6 2 )(7rV - 2tt 2 6 2 z + (1 + 7T 2 e 2 )^ 2 ) 

h ~ cs ~ n(l + z)e ' {bb) 

where we have put z = exp(— l/2e). Expanding this expression in z, we obtain the 
following expression: 

hcs = h ca -2h ca z + 0(z 2 ), (57) 

from which it is clear that for small e (and hence for very small z) h cs < h ca . For 
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example, if e = 0.05, h cs = 1.048091900 and therefore h ca - h cs = O(10" 4 ). 

When h = h cs , the zero mean solution becomes stable through a super-critical 
pitchfork bifurcation. It emits a couple of unstable solutions which disappear in a 
blue-sky catastrophe, i.e. in a collision with the stable solutions ([54)1 . These unstable 
solutions are sketched in Figs. El for various values of h, h cs < h < h ca . 

We summarise the discussion of the inclusion ()52j) in a theorem 

Theorem 4.2. For h satisfying h cs < h < h ca , the inclusion \5ty) has five 1-periodic 
solutions. 

4-3. The case of large (3 

It is possible to extend the existence of 5 solutions from j3 = oo case to the case of large 
(3 by continuity. In our numerical experiments we found that (3 does not have to be very 
large. Figure contains the plot of 3 stable solutions of equation (JU) with (3 — 3, e = 0.05, 
h = 0.544 and initial conditions x(0) = -0.413600598, -0.40456237, 1.491606074. 

Theorem 4.3. For any h satisfying h cs < h < h ca and sufficiently large (depending on 
h) f3 equation has at least five periodic solutions. 

The proof of this theorem is based on an upper semicontinuity result for differential 
inclusions j2]. To state it, we need to introduce some notation. Let B be the unit ball 
in W 1 with centre at the origin. Consider a differential inclusion 

x'(t) G F(t, x(t), A), (58) 

{t,x,X) G [0, T] x Q x A, where Q is a bounded subset of R n , A C M and F is an 
upper semicontinuous set- valued map with uniformly bounded compact images. Assume 
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that for all A G A and any x in some set Q G H the solutions of the differential 
inclusion with x(0) = Xq exist on the interval [0, T] and remain in Q. The Poincare map 
Vl{x Q ): Q -> W 1 is the set- valued map defined by 

P^(x ) = {x{T):x'{t) G F(t,x(t), A),x(0) = x }. (59) 

The following theorem is a slight variation of the continuity theorems in Section 2.2 of 
® 

Theorem 4.4 (Dependence on Parameter). Let the map (t,x,X) — > F(t,x,X) 
satisfy the conditions of the previous paragraph. Then for any Ao G A and any open 
set U C M, 2n , G U , there is a 5 > such that |A — A | < 6 implies 

graph(Vl) C graph{Vl ) + [/. (60) 

Proof of Theorem Clearly, F(t,x,/3) defined by the right-hand side viewed as 
a differential inclusion satisfies the conditions of Theorem 14.41 since all solutions of this 
inclusion are globally defined and remain bounded. 

Fix h G (h cs , h ca ). By theorem 14.41 f° r sufficiently large f3 the graph of the Poincare 
map Vp of the Suzuki-Kubo equation (jlj) is contained in a neighbourhood of the graph 
of Poo. By Theorem 14.21 the graph has five intersections with the diagonal. Therefore, 
for sufficiently small neighbourhood, i.e. for sufficiently large /3, we get at least five 
1-periodic solutions of (jlj). □ 

5. Conclusions 

In this paper we have resolved a minor controversy in the physical literature. However, 
some major questions still remain. One we can formulate as a conjecture: 

Conjecture 1. An ODE of the form of (0) with f(y) bistable has at most five 1-periodic 
solutions. 

The other fundamental question is to formulate necessary conditions on f(y) for 
(JBJ) to have a specified maximum number of solutions. We end with two remarks. 

(i) Consider (JHJ in which the nonlinearity f(y) is piecewise linear, 

(-l-y if y< -1/2 
f(y)=<V if y G [-1/2, 1/2] (61) 

(l-y if y > 1/2. 

Then it can be shown, by explicitly constructing the stable periodic solutions, that 
there is a range of h for which (jUJ) with such a nonlinearity will have five 1-periodic 
solutions jlj. This perhaps indicates that there may be more than one mechanism 
for the appearance and disappearance of solutions. 

(ii) A possible way to attack the problem of necessary conditions is to consider a 
homotopy of (J2J) into (JTJ), i.e. 

x' = (1 - 2a)x + (1 - a)tanh(/3x) - ax 3 + /icos(2vrt), (62) 
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and see, at least numerically, for which value a one no longer has a subcritical 
bifurcation from the symmetric solution. 
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Appendix A. Calculating h cs 

The definition of h cs is the value of h for which there exists a periodic solution x c (t) of 
equation 

dx 

e— = 1 - x + /icos(2vrf), (A.l) 

with the property that it connects the points (0, —1) and (0, 1) (see diagram). In other 
words, the following equations are satisfied 

/icos(27rt ) =1, (A.2) 
x c {t ) = 0, (A.3) 

x c (t + 1/2) = (A.4) 

The general solution of equation (|A.1|) is 

(A-n -*A i 1 + 4?r2e2 + h CQS ( 27rt ) + 2h7re sin (27rt) 
x[t) - Oie + - - [A.b) 



From (lA~2l we get hsm(2%t ) = y/h 2 - 1. Equations ([A~3j) - (fA~4jl then read 



1 + 4iT z e z 



Cie ^ + «v^EI = (A . 7) 

We eliminate Cie - '^ to get 

e -i/ 2e ^ 2 + 47r 2 e 2 + 27 rev / ^ r l) - 47r 2 e 2 + 27re V / ^ 3 T = (A.8) 

Denoting e _1,/2e = we finally obtain 

y/(l + 47r 2 e 2 )(7r 2 e 2 - 27T 2 e 2 ^ + (1 + 7r 2 e 2 )^ 2 ) 



7re(l + 2) 



(A.9) 



References 



[1] M. Acharyya and B. K. Chakrabarti, Magnetic hysteresis loops as Lissajous plots of relaxationally 
delayed response to periodic field variation, Physica A 202 (1994), 467-481. 

[2] J. -P. Aubin and A. Cellina, Differential Inclusions, Springer- Verlag, Berlin 1984. 

[3] N. Berglund, Adiabatic Dynamical Systems and Hysteresis, Ph. D. thesis, Institut de Physique 
Theorique, EPFL, Lausanne 1998. 



Multiplicity of periodic solutions in bistable equations 



18 



[4] G. Berkolaiko, Lecture at University of Strathclyde Analysis Seminar, 1997. 
[5] J. Byatt-Smith, Unpublished notes, 2002. 

[6] M. Golubitsky and D. G. Schaeffer, Singularities and Groups in Bifurcation Theory, Vol.1, 

Springer- Verlag, New York 1985. 
[7] A. Hohl et al, Scaling laws for dynamical hysteresis in a multidimensional laser system, 

Phys. Rev. Lett. 74 (1995), 2220-2223. 
[8] P. Jung, G. Gray, R. Roy, and P. Mandel, Scaling law for dynamical hysteresis, Phys. Rev. Lett. 

65 (1991), 1873-1876. 

[9] G. Korniss, P. A. Rikvold and M. A. Novotny, Absence of first-order transition and tricritical 

point in the dynamic phase diagram of a spatially extended bistable system in an oscillating 

field, Phys. Rev. E 66 (2002), 056127. 
[10] V. A. Pliss, Nonlocal Problems of the Theory of Oscillations, Academic Press, New York 1966. 
[11] M. Rao, H. R. Krishnamurthy, and R. Pandit, Magnetic hysteresis in two model spin systems, 

Phys. Rev. B 42 (1990), 856-884. 
[12] M. Suzuki and R. Kubo, Dynamics of the Ising model near the critical point. I, J. Phys. Soc. Japan 

24 (1968), 51-60. 

[13] T. Tome and M. J. de Oliveira, Dynamic phase transition in the kinetic Ising model under a 

time-dependent oscillating field, Phys. Rev. A 41 (1990), 4251-4254. 
[14] M. F. Zimmcr, Ising-model in an oscillating magnetic field: mean-field theory, Phys. Rev. E 47 

(1993), 3950-3955. 



